Kinetic description of particle emission from expanding source 
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The freeze out of the expanding systems, created in relativistic heavy ion collisions, is discussed. 
We combine kinetic freeze out equations with Bjorken type system expansion into a unified model. 
The important feature of the proposed scenario is that physical freeze out is completely finished in 
a finite time, which can be varied from (freeze out hypersurface) to oo. The dependence of the 
post freeze out distribution function on the freeze out time will be studied. As an example, model is 
completely solved and analyzed for the gas of pions. We shall see that the basic freeze out features, 
pointed out in the earlier works, are not smeared out by the expansion of the system. The entropy 
evolution in such a scenario is also studied. 

PACS numbers: 24.10.Nz, 25.75.-q 



I. INTRODUCTION 



In the ultrarelativistic heavy ion colhsions at RHIC the 
total number of the produced particles exceeds several 
thousands, therefore one can expect that the produced 
system behaves as a "matter" and generates collective 
effects. Indeed strong collective flow patterns have been 
measured at RHIC, which suggests that the hydrody- 
namical models are well justified during the intermediate 
stages of the reaction; from the time when local equi- 
librium is reached until the freeze out (FO), when the 
hydrodynamical description breaks down. During this 
FO stage, the matter becomes so dilute and cold that 
particles stop interacting and stream towards the detec- 
tors freely, their momentum distribution freezes out. The 
FO stage is essentially the last part of a collision process 
and the main source for observables. 

In simulations FO is usually described in two extreme 
ways: either FO on a hypersurface with zero thickness, or 
FO described by volume emission model or hadron cas- 
cade, which require an infinite time and space for a com- 
plete FO. At first glance it seems that one can avoid trou- 
bles with FO modeling using hydro-|-cascade two module 
model [H, since in hadron cascades gradual FO is real- 
ized automatically. However, in such a scenario there 
is an uncertain point, actually uncertain hypersurface, 
where one switches from hydrodynamical to kinetic mod- 
eling. First of all it is not clear how to determine such 
a hypersurface. This hypersurface in general may have 
both time-like and space-like parts. Mathematically this 
problem is very similar to hydro to post-FO transition 
on the infinitely narrow FO hypersurface, therefore all 
the problems discussed for FO on the hypersurface with 
space-like normal vectors will take place here. Another 
complication is that while for the post FO domain we 
have mixture of non-interacting ideal gases, now for the 



hadron cascade we should generate distributions for the 
interacting hadronic gas of all possible species, as a start- 
ing point for the further cascade evolution. The volume 
emission models are also based on the kinetic equations 
d, H, 0, Q defining the evolution of the distribution func- 
tions, and therefore these also require to generate initial 
distribution functions for the interacting hadronic species 
on some hypersurface. 

In the recent works Q it has been shown that the 
basic assumptions of the Boltzman Transport Equation 
(BTE) are not satisfied during FO and therefore its de- 
scription has to be based on Modified BTE, suggested 
in [6]. When the characteristic length scale, describing 
the change of the distribution function, becomes smaller 
than mean free path (this always happens at late stages 
of the FO), then the basic assumptions of the BTE get 
violated, and the expression for the collision integral has 
to be modified to follow the trajectories of colliding par- 
ticles, what makes calculations much more complicated. 
In fact, in cascade model one follows the trajectories of 
the colliding particles, therefore what is effectively solved 
is not BTE, but Modified BTE. Once the necessity of the 
Modified BTE and at the same time the difhculty of its 
direct solution were realized, the simplified kinetic FO 
models become even more important for the understand- 
ing of the principal features of this phenomenon. Such 
simplified models were developed and studied in Refs. 
[E 0) B B but all these were missing expansion. 
Thus, the important question to be studied is whether 
the important freeze out features, pointed out in the ear- 
lier works, are not smeared out by the expansion of the 
system? 

In this paper we present a simple kinetic FO model, 
which describes the freeze out of particles from a Bjorken 
expanding fireball [ll[. Thus, this is a more physical ex- 
tension of the oversimplified FO models without expan- 
sion d, 0, H, H, • Taking the basic ingredients of the 
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FO simulations from Refs. [9|, |10| we can make physi- 
cal freeze out to be completely finished in a finite time, 
which can be varied from (freeze out hypersurface) to 

00. In the other words our freeze out happens in a layer, 

1. e. in a domain restricted by two parallel hypersurfaces 
T = Ti and T = Ti + L (r is the proper time). 

In Ref. [8|] authors have also adopted kinetic gradual 
FO model to Bjorken geometry, but combined it with 
Bjorken expansion on the consequent, not on the par- 
allel basis: system expands according to Bjorken hydro 
scenario, but when it reaches beginning of the FO pro- 
cess system stops expansion and gradually freezes out in 
a fixed volume. However, it was shown in Q that such 
a scenario is not physical: the simultaneous modeling of 
expansion and freeze out is required in order to avoid 
decreasing total entropy. Now we propose such a gener- 
alized model. 



II. FINITE LAYER FREEZE OUT 
DESCRIPTION 

Let us briefly review gradual FO model, which we are 
going to generalize including expansion. Many building 
blocks of the model are Lorentz invariant and can be ap- 
plied to both time-like and space-like FO layers, so at 
the beginning we will write these in general way. Start- 
ing from the Boltzmann Transport Equation, introducing 
two components of the distribution function, /: the in- 
teracting, and the frozen out, ones, (/ = /* + f^), 
and assuming that FO is a directed process (i.e. ne- 
glecting the gradients of the distribution functions in the 
directions perpendicular to the FO direction comparing 
to that in the FO direction) we can obtain the following 
system of the equations [1, Q : 

Pesc ji ^ feq{s) — f df ^ Pesc j-i ^^-^ 

ds Tpo Tth ' ds Tpo 

The FO direction is defined by the unit vector dcr^. FO 
happens in a layer of given thickness L with two par- 
allel boundary hypersurfaces perpendicular to dcr^^, and 
s — da^^x'^'' is a length component in the FO direction. 
We work in the reference frame of the front, where da^ 
is either (1, 0, 0, 0) for the time-like FO, or (0, 1, 0, 0) for 
the space- like FO. The tfo is some characteristic length 
scale, like mean free path or mean collision time for time- 
like FO. The rethermalization of the interacting compo- 
nent is taken into account via the relaxation time approx- 
imation, where fi approaches the equilibrated Jiittner 
distribution, /eg(s), with a relaxation length/time, Tth- 
The system ((T|) can be solved semi-analytically in the 
fast rethermalization limit [oi]. 

The basis of the model, i.e. the invariant escape rate 
of the particles within the FO layer of the thickness L, 
for both time-like and space-like normal vectors is given 
as (see Refs. d, [Tol . [T2I for more details) 



where is the particle four-momentum, is the flow 
velocity. In fact the model based on the escape rate ^ 
is a generalization of a simple kinetic models studied in 
Refs. [i,0,[i, which can be restored in the L — > 00 limit. 
Here we will concentrate on the time- like case only, where 
the above O function is unity. 

The important feature of the escape rate in the form 
([2) is that physical freeze out is completely finished when 
s — i.e. it requires finite space/time extent. Further- 
more, now we can vary this layer thickness, L, from 
(freeze out hypersurface) to 00 and study how the post 
FO distribution depends on the layer thickness. Interest- 
ing and unexpected result was found in 0, [l^, [l^], for 
both space-like and time-like FO layers, namely that if L 
is large enough, at least several Tpo, then post FO dis- 
tribution gets some universal form, independent on the 
layer thickness. 

Simple semianalytically solvable FO models studied in 
@i .7; 8, j3, .IQ] are missing an important ingredient - the 
expansion of the freezing out system. The open question 
is whether the features of the FO discussed above will 
survive if the system expansion is included. In this work 
we are going to build a model, which includes both grad- 
ual FO and Bjorken-like expansion of the system, and 
answer this question. 

III. BJORKEN EXPANSION WITH GRADUAL 
FREEZE OUT 

Let us first remind the reader the basics of the Bjorken 
model Bjorken model is one-dimensional in the same 
sense as discussed before eq. ([T]) - only the proper time, 
T — — x^, gradients are considered. Here the refer- 
ence frame of the front, where da^ = (1,0,0,0), is the 
same as the local rest frame (comoving frame), where 
= (1, 0, 0, 0). The evolution of the energy density and 
baryon density is given by the following equations: 

de e + P dn n 
dT T ' dT r ' 

where P is the pressure. The initial conditions are given 
at some t — tq: e(ro) = gq, n(To) = tiq. This system can 
be easily solved: 

e(r) = eo(^)''"° , n(r) = no (^) , (4) 

where P — c^e is the equation of state (EoS) with con- 
stant velocity of sound, Cq. 

It is important to remember that if we want to have a 
finite volume fireball, we need to put some boarders on 
the system. Here we assume that our system, described 
by the Bjorken model, is situated in the spacial domain 
1*7 1 < VR or what is the same \z\ < zrIt) = Tsiuinjn 
(?7 = ^ In (jz^^ is a space-time rapidity). Within these 
boarders system is uniform along r = const hyperbolas 
due to model assumptions, while outside we have vacuum 
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with zero energy and baryon densities as well as pressure. 
Thus, we have a jump, a discontinuity on the boarder, 
which stays there during all the evolution. Certainly, to 
prevent matter expansion through such a boarder (due 
to strong pressure gradient) some work is done on the 
boarder surface [l3l |. One can think about it as putting 
some pressure to the surface with the vacuum, exactly 
the one which would remove discontinuity, then work is 
done by the expanding system against this pressure. As 
the system expends the volume of the fireball increases 
as 



V{t) = 2Axy sinh rj RT , 



(5) 



where A^y is the transverse area of the system. It has 
been shown in that for such a finite Bjorken system 
the total energy is given by 



E{t) = 2AxySmhriRTe{T) 



(6) 



while the total entropy and the conserved charge appear 
to be 

Sir) = 2A^yrjRTsiT) , N{r) = 2A^yr]RTn{T) . (7) 
Work done by the expanding system, W , is given by 

(8) 



dW = PdV 



W{t) = eoV{To) 1 - 



One can easily check then the energy conservation: 

E{t) + W{t) = eoV{To) = E{tq) = const . (9) 
Applying our FO model to such a system, we obtain: 

C(r') = - — j^nr') + ^ [Ur') - fir')] , 

(10) 

(11) 



TFO L-t' 



where FO begins at r = ri and t' = r — ri. Taking 
the fast rethermalization limit, similarly to what is done 
in Q, we can obtain simplified equations for which 
is a thermal distribution /'(r) = feq{T), and /■'', and 
consequently for e*, n* and : 



dr' Tpo L — t' ' dr' 

de^ e* L dn^ 

dr' Tpo L — t' ^ dr' 



TFO L-t' 



L 



TFO L~t' 



(12) 



(13) 



Now the idea is to create a system of equations which 
would describe a fireball which simultaneously expands 
and freezes out. Let us put our two components (e — 



+ e^) into the first equation of Q and do some simple 
algebra: 



de^ deJ_ 
dT dT 



L 



L 



T T tfo L — t' tfo L — t' 

(14) 

where last two terms add up to zero; the free component, 
of course, has no pressure. So far our eq. (fT4|) is com- 
pletely identical to the first equation of ([3]). The main 
assumption of our model is that our system evolves in 
such a way that eq. pi)) is satisfied as a system of two 
separate equations for interacting and free components 
Il5|,a6i]: 



de^ 
d7 



L 



Tfo L + ti - t 



deJ 



eJ 



L 



T TFO L + Tl-T 



(15) 



(16) 



Similarly we can obtain equations for baryon density [ISl . 



dn^ 
dn^ 



T Tfo L + Ti — t 



T 



L 



Tfo L + Ti-T 



(17) 



(18) 



In all these eqs. (|15m8|) on r.h.s. there are two terms: 
the first one is due to expansion and the second one is 
due to freeze out. 

As it was discussed above our volume is restricted in 
77 space by two sharp boarders, |t7| < r]R = const. All 
our interacting matter is maintained within these limits. 
However, the frozen out particles, which do not interact 
with system anymore, will not necessarily know about 
such a restriction and can, in principle, move faster than 
the boarder and thus will leave the volume under consid- 
eration. The momenta of the frozen out particles are dis- 
tributed according to the distribution function, /■'', which 
(as we shall see later) falls down exponentially with mo- 
mentum, while only the particles strongly directed along 
the beam line can have high rapidity. So, from the par- 
ticles emitted close to midrapidity only a very tiny frac- 
tion will have momentum rapidity higher than rjji, and 
thus will leave our system volume. Particles freezing out 
from the system elements with higher space-time rapid- 
ity will have better, but still small, chance to leave the 
volume of consideration, and only the frozen out particles 
from the volume elements close to the \ri\ = rjR boarders 
will disappear from our modeling with probability up to 
50%, because is spherically symmetric in the comov- 
ing frame. For the reasonable r]R boarders the amount of 
such particles can be neglected. It is, however, clear that 
our approach is best justified for the reaction region near 
midrapidity. In this work we will do only illustrative cal- 
culations, aiming for a qualitative understanding of the 
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general features of the FO and checking the relation to 
the former models of this process [1, 0, [13], while for 
the more quantitative calculations and comparison with 
the experimental data we will restrict our modeling to 
midrapidity region p^ . 

Thus, finally, we have the following simple model of 
fireball created in relativistic heavy ion collision. 
Initial state, t = tq bq, uq 

Phase I, Pure Bjorken hydrodynamics, tq < t < ti 



e(r) = eo (^) 



n{r)^no{^] (19) 



Phase II, Bjorken expansion and gradual FO, 

Ti < T < Ti + L 

Solving Eqs. (|15|17p we obtain: 



eHr) = eo 



To\l+=o fL + Tl-T 



L 



nXr)=no(^) 



L + Tl 



l/tfo 



L/tfo 



(20) 



(21) 



The difference with respect to the pure Bjorken solution 
(|19p is in the last multiplier, and we see that, as expected, 
the interacting component completely disappears when r 
reaches r = L + Ti . 

With these last equations we have completely deter- 
mined evolution of the interacting component pisl . [l^. 
Knowing e*(T) and EoS we can find temperature, Ti(T). 
Due to symmetry of the system u^{t) = u^(ro) = 
(1,0,0,0). Finally, /'(r) is a thermal distribution with 
given T,(r), n*(r), <(t). 

However, for us the more interesting is the free compo- 
nent, which is the source of the observables. Eqs. (|16ll8p 
give us the evolution of the e / and n f , and one can eas- 
ily check that these two equations are equivalent to the 
following equation on the distribution function: 



dT 



Tpo L + Tl-T 



(22) 



The measured post FO spectrum is given by the distri- 
bution function at the outer edge of the FO layer, i.e. by 

//(L-Hn). 

Most of the observables will depend only on this mo- 
mentum distribution. However the two particle correla- 
tions depend also on the space-time origins of the de- 
tected particles. Thus, in order to calculate two parti- 
cle correlations we have to keep track on both the mo- 
mentum and the space-time coordinates of the freeze out 
point of particles (this can be done, for example, in the 
source function formalism), i.e. we can use the full infor- 
mation regarding the space-time evolution of free com- 
ponent, eq. (j22p . and correspondingly can get some re- 
strictions on it from the data. 



IV. RESULTS FROM THE MODEL 

Aiming for a qualitative illustration of the FO process 
we show below the results for the ideal massive pion gas 
with Jiittner equilibrated distribution T?]: 



r(T,p) = 



(27r)3 



(27r)3 



,-VlpP+™jm(r) ^ 

(23) 

in the comoving frame, where it^ — (1,0,0,0). Here the 
degeneracy of pion is g — 3, while the baryon chemical 
potential in the case of pions is zero. 

Contrary to the illustrative example in [l5l | here we do 
not neglect the pion mass. During FO the temperature 
of the interacting component decreases to zero, so at late 
stages of the FO process this new calculation is better 
justified. We will see below that Ti falls below quite 
soon, and so the Jiittner distribution is a good approxi- 
mation of the proper Bose pion distribution [Tg} . 

For our system we have the following EoS: 



27r2 



(24) 



where if„ is Bessel function of the second kind, and a — 
m/Ti. 

Eqs. and ([24]) result into the following equation 

for the evolution of the temperature of the interacting 
component: 

dT, ^ T, 4TfK2{a)+mT,Ki{a) 

dT ^ T l2TfK2{a) + 5mT,Ki{a) + w?Ko{a) 



T- 
tfo 



L 



3T^K2{a)+mT,Ki{a) 



L + Ti-T J l2TfK2{a) + ^mT,Ki{a) + m^Koia) 

(25) 

Furthermore, we have used the following values of the pa- 
rameters: rjn = 4.38, Axy = TrR\^, where Rau = 7.685 
fm is the Au radius, tq = 0.05 /m, T^{to) = 820 MeV, 
ri = 5 fm, what leads to T,{ti) = Tpo = 180 MeV, and 
Tfo = 0.5 fm. During the pure Bjorken case the evo- 
lution of the temperature is govern by eq. ([^5]) without 
the second (freeze out) term on the r.h.s. 

In Fig. [l]we present the evolution of the temperature 
of the interacting matter, T'i(T), for different values of FO 
time L, and Figs. [5] present the evolution of the energy 
densities for the interacting and free components. These 
satisfy the energy conservation Ef{T) + Ei{T) -\- W{t) = 
E{to), where Ef{T) and Ei{T) are given by eq. ^ with 
corresponding energy density. 

As it was already shown in Ref . [1, [13] , the final post 
FO particle distributions, shown on Fig. [3l are non- 
equilibrated distributions, which deviate from thermal 
ones particularly in the low momentum region. By in- 
troducing and varying the thickness of the FO layer, L, 
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T, = 5.0 fm/c, =0.5 fm/c 




T[fm] 

FIG. 1: Evolution of temperature of the interacting matter for different FO layers. Ti{To — 0.05 fm) = 820 AleV, Tfo = 
180 MeV . "No Freeze Out" means that we used standard Bjorken hydrodynamics even in phase II. 



we are strongly affecting the evolution of the interacting 
component, see Fig. [U but the final post FO distribu- 
tion shows strong universality: for the FO layers with 
a thickness of several Tpo, the post FO distribution al- 
ready looks very close to that for an infinitely long FO 
calculations, see Fig. |3]left plot. Differences can be ob- 
served only for the very small momenta, as shown in Fig. 
[3] right plot. So, the inclusion of the expansion into our 
consideration does not smear out this very important fea- 
ture of the gradual FO. 

Please note, that if one would look on our post FO 
distributions only in the medium momenta regions (with 
0.5 GeV < \p\ < 2.5 GeV), then one could fit these spec- 
tra reasonable well with equilibrated distribution with 
temperature T ~ 160 MeV, see dashed line on Fig. [31 
However, for low and high momenta such a fit would 
strongly disagree, this has been shown already in earlier 
FO model in Ref. [7b]. 

For our illustrative calculation we assumed rjn = 4.38 
what is one unit smaller than the original rapidity of the 
colliding nuclei for 100 GeV * A. The experimental data 
show a smaller plato in rapidity, but we used a classic one- 
dimensional Bjorken model without transverse expansion 
which correspondingly requires a larger longitudinal size. 

It is important to check the non-decreasing entropy 
condition [l5l . [TgI . ITsj to see whether such a process is 
physically possible. Figs. [J] present the evolution of the 
total entropy, S{t), calculated based on the full distribu- 



tion function, /(p) = /'(p) + /^(p): 



1 - In 



(2^ 
9 



fir] 



(26) 



The total entropy is then given by eq. ([7|. 

From this figure we can make an important conclu- 
sion, that gradual freeze out with rethermalization pro- 
duces entropy. This is not new - in Ref. [3] the au- 
thors reached the same conclusion based on longitudially 
boost invariant hydrodynamical model with continuous 
emission (what means infinitely long FO in our language) 
with different escape probability [2]. In Ref. [2, 3] Pesc is 
given by the probability to escape without further colli- 
sions, calculated based on Glauber formula. In the case of 
complete rethermalization of the interacting component, 
which we also assume, their entropy increase can go up 
to 70%! This leads then to the corresponding increase in 
the pion production, which is a rough measure of entropy. 
Certanly such a hydro+continuous emission model leads 
the nonthermal final distributions, and authors found a 
reasonable agreement with NA35 transverse mass spec- 
tra. This model naturally explains the large number of 
the produced pions, which was a trouble for models with 
sharp FO hypersurface. 

In our model, as well as in 0, , there are two main 
sources of entropy production, namely the separation 
process (into interacting or free component) and the 
rethermalization process. Separation process leads to en- 
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FIG. 2; Evolution of the energy densities for the interacting 
(a) and free (b) components for different FO layers during 
phase II. Ti{To = 0.05 fm) = 820 MeV, Tfo = 180 MeV. 



tropy increase because particles have now more accessible 
states being either interacting or free, not just interact- 
ing. The entropy increase in the rethermalization process 
is a manifistation of the H theorem. 

At the late stages of the long FO this entropy increase 
can be approximated analytically in our model. 



dr J dr 



1 - In 



(27r) 



we obtain 



dl 
dr 



dr 



(28) 



Then, putting eq. ([^5]) into eq. (P7)) we get the following 
equation for the entropy evolution: 



ds{T) s{t) ^ 1 
dr T 



- I d'pf 



f 



(29) 



If there is a baryon charge (or other conserved charges) 
in the system in general case the last term in the eq. ([29)) 
can then be related to rif: 



ds{T) s(r) nf (t) 



dr T 
and for the baryon density, eq. 

dn^ [t) (t) 



(30) 



TSl). in this limit we have 



nf{T) = npo- 



dT 



Tl 



(31) 



where npo is the free baryon density at the end of FO. 
So, for the total entropy, eq. ([7]), we obtain: 



dS d (st) 



S 

T 



S 

T 



N 



FO 



(32) 



where Nf{T) is the total number of frozen out baryon 
charges, which is in our limit, eq. PT|) . equal to the 
total number of frozen out baryon charges at the end of 
FO, Npo, which is, due to conservation law, the same as 
initial number of interacting baryon charges. 



S{t) = Sfo + Nfo In 



L + Tl 



(33) 



where Spo is the final entropy. 

Thus, we see that the total entropy increases during 
simultaneous expansion and gradual FO, as it is clearly 
seen on Fig. |4]for the L — > oo curve. If there is a baryon 
charge in the system, then at the very end of the FO 
process entropy increases logarithmically, according to 
eq. (IMl). 

On the other hand the fast FO in the narrow layer, 
which can be studied in our model, contrary to 0, 
does not lead to entropy increase. Already for the FO 
layer as thin as 2tfo the entropy production is negligible, 
around 1%, see Fig. ID 



V. CONCLUSIONS 



- / d^P 



df{r) 
dr 



(27) 



At the late stages of the FO, when r is close to the end 
of FO L + Ti, for long FO, i.e. L » Tpo, the interacting 
component, fi, can be neglected next to free component, 
//, (see eqs. (|20I21|) and Fig. O), and so, using eq. 



In this paper we presented a model, which describes 
simultaneously freeze out and Bjorken expansion, and 
thus, it is more a physical extension of the oversimplified 
FO models without expansion 0, H, HI, > which al- 
lowed us to study FO in a layer of any thickness, L, from 
to oo. Another important feature of the proposed model 
is that it connects the pre FO hydrodynamical quanti- 
ties, like energy density, e, and baryon density, n, with 
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p[GeV] p[GeV] 

FIG. 3: Final post FO distribution for different FO layers as a function of the momentum in the FO direction, p = in our 
case (p" = p^ = 0). The initial conditions are specified in the text. Dashed curve shows thermal distribution with temperature 
T = 160 MeV . Subplot (b) is a magnified view on the low momentum region of subplot (a). 
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FIG. 4: Evolution of the total entropy for different FO layers. The initial conditions are specified in the text. Subplot (b) is a 
magnified view on the beginning of FO region of subplot (a). 



the post FO distribution function in a relatively simple 
way. As a zero order approach to the real physical sys- 
tem we studied the FO of the pion gas [lB| . 

The results show that the inclusion of the expansion 
into FO model, although strongly affects the evolution 
of the interacting component, does not smear out the 
universality of the final post FO distribution, observed 
already in Refs. 



Another important conclusion of this work is stress- 
ing once again the importance to al way s check the non- 
decreasing entropy condition since long grad- 
ual freeze out may produce substantial amount of en- 
tropy, as shown in Fig. |4l 

This may have important consequence for QGP search. 
For many qualitative estimations it was assumed that 
all the entropy is produced at the early stages of the 
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reaction and that the expansion, hadronization and FO 
go adiabatically, and thus number of pions can serve as 
a rough measure of entropy. However, if a non-neghgible 
part of the entropy, say 10%, is produced during FO, then 
some estimations, for example of strangeness vs entropy 
(pion) production have to be reviewed. 
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